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Abstract 



Soft-collinear effective theory is used to prove factorization of the B — > 'jli' decay 
amplitude at leading power in A/mt, including a demonstration of the absence 
of non-valence Fock states and of the finiteness of the convolution integral in the 
factorization formula. Large logarithms entering the hard-scattering kernel are 
resummed by performing a two-step perturbative matching onto the low-energy 
effective theory, and by solving evolution equations derived from the renormal- 
ization properties of the leading-order i?-meson light-cone distribution amplitude. 
As a byproduct, the evolution equation for heavy-collinear current operators in 
soft-collinear effective theory is derived. 



1 Introduction 



The proposal of a "soft-coUinear effective theory" (SCET) for the strong interactions 
of colhnear and soft particles has been an important step toward understanding the 
factorization properties of hard exclusive processes in QCD [H |2l El EI- In particular, 
it raises the prospects for rigorously proving QCD factorization theorems for hadronic 
and radiative i?-meson decays into light particles such as -B ^ tttt jH] and B K*^ 
ini m, which are of great importance to the physics program at the B factories. A 
challenge common to these decays and many others is to understand the interactions of 
colhnear particles with the soft spectator quark inside the B meson, which give rise to 
convolutions of hard-scattering kernels with 5-meson light-cone distribution amplitudes 
(LCDAs). The first systematic analysis of these interactions in the framework of SCET 
has recently been performed by two of us 

The radiative, semileptonic decay B 7/1/ provides a clean environment for the 
study of soft-collinear interactions |9j. This process is particularly simple in that no 
hadrons appear in the final state. Yet, there is sensitivity to the light-cone structure of 
the B meson, probed by the coupling of the high-energy photon to the soft spectator 
quark inside the heavy meson. In the present paper, we apply the formalism of [S] to 
prove factorization for this decay and systematically resum large Sudakov logarithms. 
The arguments we will present apply, with minimal modifications, to more complicated 
decays such as i? — > K*'y. 

Several other groups have recently studied the decay B •ylu. In (TU] a QCD 
factorization formula was established at next-to-leading order (NLO) in as, and it was 
demonstrated that the leading-order LCDA of the B meson is sufficient to describe the 
decay amplitude at leading power in A/ mi, (with A a typical hadronic scale), contrary to 
the findings of p. QCD factorization formulae have also been proposed for the related 
processes i? ^ 77 and i? — 7 |TT| ll2j. Arguments in favor of factorization in higher 
orders were given in jT^ using a formulation of SCET different from the one adopted 
here. 

In the present work we provide the first complete proof of factorization for a B decay 
in which hard-spectator interactions are relevant at leading power. We believe that our 
approach is simpler and more transparent than that put forward in . The discussion of 
factorization we will present is more complete in that we prove the absence of non- valence 
Fock-state contributions and the convergence of the convolution integral to all orders 
in perturbation theory. We present the first correct result for the perturbative hard- 
scattering kernel in a scheme which uses the conventional, mass-independent definition 
of the LCDA. In addition, we discuss in detail the complete renormalization-group (RG) 
resummation of large logarithms. To this end, we solve evolution equations for the 
different components of the hard-scattering kernel, which follow from the renormalization 
properties of the B-meson LCDA. In this context we clarify the connection between 
the anomalous dimension of heavy-collinear currents in SCET and the cusp anomalous 
dimension encountered in the study of Wilson loops with light-like segments |14[ I15j. 
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Our main goal is to establish the QCD factorization formula jTU] 

A{B' ^ jrui) (X rriBfBQu [ dl+ '^+^|+'^-^ T{1+, E^,mb, fj,) (1) 

Jo 

to all orders in perturbation theory and at leading power in A/m^. Here = | is 
the electric charge of the up-quark (in units of e), fs is the 5-meson decay constant, 
(j)^ is a leading-order LCDA of the 5- meson, and T = 1 + 0{as) is a perturbative 
hard-scattering kernel. Factorization holds as long as the photon is energetic in the 
i?-meson rest frame, meaning that is of the order of the 6-quark mass. The physics 
underlying the factorization formula is that a high-energy photon coupling to the soft 
constituents of the B-meson produces quantum fluctuations far off their mass shell, which 
can be integrated out in a low-energy effective theory. Specifically, when the photon 
couples to the 6-quark it takes it off shell by an amount of order ml, producing a hard 
quantum fluctuation that can be treated using the methods of heavy-quark effective 
theory (HQET) ^Hl- (The resulting contribution to the amplitude is, in fact, power 
suppressed.) When the photon couples to a soft light parton inside the B meson it 
produces a "hard-collinear" mode that is off shell by an amount of order mf,A. Once 
these short-distance modes are integrated out, the decay amplitude factorizes into a soft 
component (the LCDA) and a hard-scattering kernel. 

For the analysis of the B ^Iv decay amplitude we work in the S-meson rest 
frame and choose the photon momentum along the z direction, such that g'^ = E^n^, 
where = (1, 0, 0, 1) is a light-like vector. The two transverse polarization states of 
the photon can be expressed in terms of the basis vectors = ^(0, 1, =F^, 0), which 
correspond to left- and right-circular polarization, respectively. It is convenient to define 
a second light-cone vector = (1,0,0,-1). Any 4-vector can then be expanded as 

= ^{p+n^^ + p^n^) + p^, where p^ = n ■ p and p^ = n ■ p. 

In order to prove the factorization formula (0) one needs to show that |17j : 

1. The decay amplitude can be expanded in powers of transverse momenta and, at 
leading order, can be expressed in terms of a convolution with the i?-meson LCDA 
as shown in (|T|). 

2. After subtraction of infrared contributions corresponding to the B-meson decay 
constant and LCDA, the leading contributions to the amplitude come from hard 
internal lines, i.e., the hard-scattering kernel T is free of infrared singularities to 
all orders in perturbation theory. 

3. The convolution integral of the hard-scattering kernel with the LCDA is convergent. 

4. Non-valence Fock states do not give rise to leading contributions. 

In the authors have discussed factorization for the decay B 7/1/, albeit without 
addressing the last two points in this list. (The same criticism applies to the treatment of 
factorization for the decay B Drr presented in ^H]-) The verification of the third point 
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is in essence a check that one has correctly identified the infrared degrees of freedom in 
the effective low-energy theory. If the integral in diverged for 0, this would 

imply that in addition to the hard-coUinear modes mentioned above one would need to 
consider long-distance contributions from internal momenta collinear with the photon 
momentum, which are off-shell by an amount of order A^. The fourth point addresses 
contributions from non-valence Fock states, which in the present case correspond to 
i?-meson matrix elements of trilocal qG^'^h operators. While local operators of this 
type would be of higher dimension and hence power suppressed, non-local operators can 
contribute at leading power if their component fields are smeared over domains of size 
1/^0- following section we will address all of the above and, on the example 

01 B ^ 7/z/ decay, provide the first complete proof of a QCD factorization theorem 
including hard-spectator interactions. 

Before entering the technical details of the factorization proof, we stress that (P) 
provides the simplest example of a factorization formula for a decay in which the hard- 
scattering mechanism involves three different mass scales: a hard scale ~ m;,, a soft 
scale Z+ ~ A, and an intermediate scale ~ i/mfcA. Indeed, a diagrammatic 

analysis of one-loop diagrams for the decay amplitude reveals that leading contributions 
arise from three different regions of loop momenta: hard momenta k ~ rrify, soft momenta 
/c ~ A, and hard-collinear momenta scaling like {k+, k-, k±) ~ (A, rub, y^rribA). This 
suggests a second stage of "perturbative" factorization [Hllinilini, which we will establish 
below. It says that the hard-scattering kernel itself can be factorized as 

T(/^, E,, m„ ^) = H f ^, ^) . J f ^) . (2) 

\ IJ, nib J \ fi^ J 

The hard component H accounts for the short-distance corrections from quantum fluc- 
tuations that are off shell by an amount of order m^. They arise from the coupling of 
hard or hard-collinear particles to the heavy quark. The jet function J accounts for 
short- distance fluctuations that are off shell by an amount of order rubA, which result 
from the coupling of hard-collinear particles to the light spectator quark. The factoriza- 
tion formula (j21) thus separates the physics on two different short-distance scales. While 
this is necessary to gain full control over large logarithms arising in the perturbative 
calculation of the hard-scattering kernel, it is not a necessary step in the proof of the 
QCD factorization formula (Q), which describes the separation of short- and long-distance 
physics. 

2 Proof of factorization 

We use the formulation of SCET developed in [S], where it was argued that the discussion 
of exclusive B decays into light particles can be made most transparent by matching full 
QCD amplitudes onto a low-energy effective theory in which only long-distance modes 
are kept as dynamical degrees of freedom. In the present case the B meson is the only 
hadron in the process, and so the relevant degrees of freedom in the low-energy effective 
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theory are soft partons. The strong-interaction Lagrangian consists of the ordinary QCD 
Lagrangian for hght quarks and gluons (restricted to the subspace of soft Fourier modes) 
and of the HQET Lagrangian for heavy quarks. In the full theory the hadronic part of 
the decay amplitude is given by the i?-meson matrix element of the time-ordered product 
of a weak, flavor-changing current and the electromagnetic current. At leading order in 
the effective theory this object is matched onto flavor-changing bilocal operators of the 
form A^'^^^\z) qs{z) . . . h(0), where A^'^'^^ is the photon field, h is the effective-theory field 
for the heavy 6-quark, and is the field for a soft M-quark. The separation z between the 
fields is nearly light-like, ~ l/(mfeA) 0. A space-time picture of the decay process is 
as follows: The virtual light quark produced at the weak vertex is off shell by an amount 
of order mf,A and almost coUinear with the photon direction. The electromagnetic vertex 
is thus located at a distance z'^ = tn^+sn^+z^ relative to the weak vertex, where t ~ 1/A 
is the parametrically largest component. Because of the scaling properties of the photon 
and soft spectator momenta one can replace qs{z) ~ qs{tn) and A^'^™\z) ~ A^^^\sn) 
to leading power. From the discussion in it follows that at leading order in A/rrif, 
there exists a unique type of SCET operators that mediate this decay and are allowed 
by gauge and reparameterization invariance. They are 



5^ I d'x T { [uY{l - 75)fe] (0), [g4^'"^ q] 

q=u,b 

dsdtd,{t,s,vq,mb,^i) Q,(tn) 47^^^) | J{(0) 

= J2j dtCi{t,n-q,v-q,mb,fi)qs{tn)S{tn,0)4'~^r\^)^^iK0), (3) 

where "K = S'^h and = S'^qs with the path-ordered exponential 

S{x)=Pexp^g J dw n ■ As{x + wn)^ (4) 

are gauge-invariant combinations of SCET fields and Wilson lines, and Ag is the soft 
gluon field. The combination S{x,y) = S{x)S''{y) appearing in the last line of 
represents a soft Wilson line connecting the two points x and y on a straight segment. 
Physically, this string operator arises because the virtual quark propagating between 
the two vertices can emit multiple soft gluons without power suppression. The quantity 
^[.7'' is the electromagnetic analog of the gauge-invariant coUinear gluon field defined in 
jH]. To first order in e we have 

47)(0) = n,7M rfweF"^(w;n). (5) 



oo 



The Feynman rule for this object is simply e^*, where e is the photon polarization vector. 
Finally, from the fact that the leptonic weak current 2^7/^(1 — 75) / is conserved (in the 
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limit where the lepton mass is neglected) it follows that the relevant Dirac structures Fj 
in can be taken as Fi = 7'^(1 — 75) and F2 = n^{l + 75). 

In the last step in (jHl) we have used that the photon momentum q is an external 
momentum, so that the integration over s can be performed and leads to new coefficients 
Ci{t,n ■ q,v ■ q,mb,fi) = J ds e^^'^''^ Ci{t, s,v ■ g, mb,/i). Note that the scalar products 
n ■ q = 2E^ and v ■ q = are both determined in terms of the photon energy; however, 
for a while it will be useful to distinguish between these two variables. 

SCET power counting shows that the soft fields IK and Q^^cale like A^/^, the inte- 
gration variable t scales like 1/A, and the Wilson coefficients Ci scale like 1. It follows 
that the hadronic components of the SCET operators on the right-hand side of (jH)) scale 
like A^, which is one power of A less than a local current containing two soft quark fields. 
That one can write down such operators is a consequence of the fact that transverse 
collinear fields have unsuppressed interactions with soft light quarks |H|. As a result, 
at leading power the Wilson coefficients Ci receive contributions only from Feynman 
diagrams with a photon attached to the light spectator quark in the B meson. 

Following [20J, we define the two leading-order LCDAs for the B meson in position 
space in terms of the HQET matrix element 



^ '0\qs{z)S{z,0)Th{0)\B{v)) 



tr 



0?(r,/i)-^ 0^(r,/x)-0^(r,/i) 



(6) 



Here = is a null vector, v is the i?-meson velocity, and r = v ■ z — iO. The quantity 
F{fi) corresponds to the asymptotic value of the product fBy/f^B in the heavy-quark 
limit ,2.1j. The two distribution amplitudes obey the normalization 0^(0, /i) = 1 at r = 0. 
In our case z = tn, and because of the presence of the factor ^ in (jS)) it follows that only 
the function (j)^{T,fi) contributes, to which we refer as the leading-order LCDA. Note 
that with our choice of the light-cone basis vector n we have t = t v ■ n = t. Performing 
the relevant traces over Dirac matrices, we find that at leading power in A/mb the decay 
amplitude vanishes if the photon has right-circular polarization, while for a photon with 
left-circular polarization it is given by 

V2 



X ^/m^F{^l) J citCi(t,n ■ g, w ■ g, mb,/i) + . . . , (7) 

where the dots represent power-suppressed contributions. Only the SCET operator with 
Dirac structure Fi = 7^(1 — 75) contributes to the decay amplitude. In order to recast 
the above result in a form resembling the factorization formula (0) we introduce the 
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Fourier transforms of the Wilson coefficient function and the LCDA as jHl 1201 



Ci{l+,n- q,v ■ q,mb,fi) = dte n ■ g, w ■ g, m,,, /x) , 

(8) 

The analytic properties of the function /i) in the complex t-plane imply that 

(j)^{u!,fi) = if cij < 0. Finally, the HQET parameter F{^) is related to the physi- 
cal _B-meson decay constant through fsy/fnB = Kp{mi), ji) F^jj) [1 + ©(A/mf,)], where 
at NLO in the MS scheme j21j 

AV(m..^) = l + ^!:^(31n!=i-2). (9) 

Combining these results, it follows that the terms shown in the second line of (|7j) equal 
those on the right-hand side of the factorization formula ^ if we identify the hard- 
scattering kernel as 

^T{l+,E^,mb,iJ.) = Kp^imb.jj) Ci{l+,2E^, E^,mb, fi) . (10) 

Let us now discuss factorization in the context of our formalism. The fact that the 
position-space SCET operators appearing on the right-hand side of (jSl) contain compo- 
nent fields with light-like separation implies that transverse parton momenta can be set 
to zero at leading power. As we have seen, this naturally leads to the appearance of 
LCDAs. In SCET the hard-scattering kernel T is identified with a Wilson coefficient. 
The absence of infrared singularities then follows from the very existence of a low-energy 
effective theory, because Wilson coefficients arise from matching and by construction are 
insensitive to infrared physics. At this point we have achieved as much as H^. 

We proceed to prove the convergence of the convolution integral in (^. The key 
ingredient here is to note that the invariance of SCET operators under reparameteriza- 
tions of the light-cone basis vectors n and n j22| can be used to deduce the dependence 
of Wilson coefficient functions on the separation t between the component fields of non- 
local operators [S]. In our case, invariance of the operators in Q under the rescaling 
transformation ra'^ —>■ n^/a and an^ (with fixed v) implies that 

Ci{t, n- q,v ■ q, rrib, fi) = Ci{at, an ■ q,v ■ q, nib, /i) 

= Hi{v ■ q, nib, /u) ■ J{at, an ■ q, /x) (11) 

to all orders in perturbation theory. In other words, the variables t and n ■ q can only 
appear in the combination n-q/t, but not individually. (Similarly, the variable t enters the 
argument of the LCDA in (jHI) in the reparameterization-invariant combination r = tv-n.) 
In the second step we have used the fact that v ■ q and enter only through interactions 
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of hard or hard-coUinear gluons with the heavy quark. The corresponding modes can be 
integrated out in a first matching step and lead to the functions if,, which depend on 
the Dirac structure of the weak current containing the heavy quark. The non-locahties 
of the component fields in Q result from the coupling of hard-collinear fields to the soft 
spectator quark in the B meson. These effects live on scales of order m^A and can be 
integrated out in a second step, leading to the function J. Since the coefficients Q are 
dimensionless, it follows that (with a slight abuse of notation) 

Ci{t, n- q,v ■ q, nib, /i) = Hi (^-j^^ ^-r^ " J {j^^ ' ^^'^'^ 

where x-y = 2v ■ q/rrib = 2E^/mb is a scaling variable of order 1. Corrections to this 
perturbative factorization formula are suppressed by a ratio of the intermediate scale 
n ■ q/t ~ nibA and a hard scale of order m^. (While the fact that these corrections scale 
like A/rub is an immediate consequence of our discussion here, it is not obvious in the 
context of the approach proposed in [Tallin], which is based on an expansion in powers of 
A/nib.) The corresponding result for the hard-scattering kernel obtained after Fourier 
transformation has the form shown in (0) if we identify 



(13) 



where n ■ ql^ = 2q ■ I = 2E^ Since the dependence of the coefficient functions on the 
renormalization scale is logarithmic, it follows that to all orders in perturbation theory 
the Wilson coefficients in (|T^ scale like Cj ~ 1 modulo logarithms. (Correspondingly, 
the kernel scales like T ~ 1 modulo logarithms.) The convergence of the convolution 
integral in ((Tj) in the infrared region t — oo, corresponding to the region — ^ in the 
factorization formula (P), then follows to all orders in perturbation theory as long as the 
integral converges at tree level. Because the B meson has a spatial size of order 1/A due 
to confinement, the bilocal matrix element must vanish faster than 1/t for t ^ 1/A, and 
so the integral over t is convergent. 

The absence of endpoint divergences in convolution integrals is tied to the question 
of whether the infrared degrees of freedom have correctly been identified in the effective 
theory ^]. Let us illustrate this point on the case at hand. Assume the integral over 
/+ in diverged logarithmically for /+ — 0. This would indicate that there were 
leading-power contributions to the decay amplitude from regions where intermediate 
propagators (such as the quark propagator in Figure Q below) are collinear with the 
external photon. In addition to the modes we call hard-collinear one then had to include 
other, long-distance collinear modes in the construction of the effective theory, which 
would invalidate factorization. The decay amplitude would then be sensitive to the 
hadronic structure of the photon as well as the full Bethe-Salpeter wave function of the 
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B meson. Our demonstration that the region /_|_ — > is power suppressed is thus an 
essential ingredient of the factorization proof. 

The final step in the factorization proof is to demonstrate the power suppression of 
more complicated projections onto the B meson involving higher Fock states or transverse 
parton momenta. The fact that such projections do not contribute at leading power 
follows from the rules for constructing SCET operators out of gauge-invariant building 
blocks, as explained in jH]. Projections sensitive to transverse momentum components 
contain extra derivatives and so are power suppressed. Projections corresponding to 
non- valence Fock states contain insertions of the soft (subscript "s" ) gluon field 

yi^(x) = [S^{iD^S)] (x) = / dwn^ [S^gC^^S] (x + wn) . (14) 

Since As scales like A, such insertions lead to power suppression unless this field is 
integrated over a domain of extension 1/A. It follows that the only possibility for a 
leading-power contribution from non-valence Fock states would be to include a factor 
of J duA'^{un) with m ~ 1/A (up to dimensionless Wilson coefficients of order 1). An 
example of such an HQET operator (omitting gauge strings and Dirac structures) is 
J dt f dw w (jsitn) UaG^'' {wn) h{0), which is indeed of leading order ~ A^ in power count- 
ing. Reparameterization invariance requires that this object must be accompanied by 
an additional factor of n in the numerator.^ Since n ■ = by definition, the only 
possibility would be to include an insertion of J dui/iJ^si.{un) somewhere between the 
light-quark field and the Dirac matrix Pj in the SCET operators in Q. However, 
any such insertion vanishes, since = 0. It is important for this argument that the 
heavy quark can be integrated out before one removes the off-shell modes resulting from 
soft-collinear interactions. The heavy quark therefore decouples from such interactions. 
This ensures that the factor i/i cannot appear to the right of Pj, and it excludes the 
appearance of f ■ n instead of i/i. 

We have thus completed the proof of the factorization formula (0) to all orders in per- 
turbation theory, and at leading power in A/m;,. The remainder of this paper is devoted 
to the calculation of the hard-scattering kernel at NLO in RG-improved perturbation 
theory, including a complete resummation of large logarithms. 

3 Calculation of the hard-scattering kernel 

The Wilson coefficients Ci in Q are derived by matching perturbative expressions for 
operator matrix elements in the full theory onto corresponding expressions in the effective 
theory. Because by construction the Wilson coefficients are insensitive to infrared physics 
the matching can be done using on-shell external quark states. We thus assign incoming 

^The reason is that the integration variable u scales like n, and that the basis vectors n and fi always 
come in pairs since they are introduced through the light-cone decomposition of 4-vectors. The only 
alternative to having a factor of n in the numerator would be to add a factor of n/{tn- q) or n/(uh- q), 
both of which scale like A/rrib and thus would lead to power suppression. 
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Figure 1: Tree- level matching calculation for the Wilson coefficients Ci and 
C2, without (left) and with (right) an external soft gluon. The dashed line 
denotes the flavor-changing weak current. The resulting non-local operators 
in SCET are denoted by a crossed circle. Double lines represent effective 
heavy-quark fields in HQET. 

momenta m^v to the heavy quark and / (with P = 0) to the soft light quark. At 
tree-level the relevant amplitude in the full theory is obtained from the first Feynman 
diagram shown in Figure Q At leading power only photon emission from the light 
spectator quark contributes; emission from the 6-quark is suppressed by one power of 
A/rrib. The corresponding amplitude in the effective theory follows from the second 
diagram. A straightforward matching calculation yields for the Wilson coefficients at 
this order Ci = Qu/l+ and C2 = 0, where /+ = n ■ I — iQ. (In general, the variable 
/+ is conjugate to the coordinate t in Q. In the present case, /+ coincides with the 
plus component of the spectator momentum because of the particular external state we 
choose for the matching calculation.) The corresponding results in position space are 
Ci = iQu 0{t) and C2 = 0. 

While these results are most easily derived by matching amplitudes with two external 
quarks, the Wilson coefficients are independent of the nature of the external states. 
Alternatively, therefore, they can be determined by matching amplitudes with external 
soft gluons. Gluon emissions from the external quark lines cancel in the matching, since 
those emissions are the same in the full theory and in SCET. However, gluon emissions 
from the internal quark propagator, which is integrated out in the effective theory, are 
contained in the Wilson line S{tn,0), which sums up an infinite number of soft gluon 
insertions. The one-gluon example is illustrated in the last two diagrams in Figure ^ 
Evaluating these graphs one readily recovers the results given above. In that way one 
con&ms (at tree level) our general result about the absence of operators containing 
additional insertions of the soft gluon field Ag. 

Beyond tree level the coefficient functions can be written in the form 



To obtain the NLO corrections Cj we evaluate the one-loop contributions to the decay 
amplitude in the full theory and in SCET. The relevant diagrams in full QCD are shown 
in Figure |21 In addition there is a contribution from the wave-function renormalization 
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(15) 
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Figure 2: 
to the B 
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One-loop diagrams in the full theory contributing at leading power 
'-flu decay amplitude. 



for the heavy quark. Using the MS regularization scheme with d = 4 — 2e dimensions 
and anti-commuting 75, we find the following expressions for the contributions Ai corre- 
sponding to the coefficients q (before subtraction of the pole terms): 



^QCD ^ , ^ 



-2e 



1 

2e 



In 



2 l-\ 



2(1 - 21ogx^ 



In^ 



4 In^ 



mt 



+ 



2 3^/^ 



Inx-, — 2^2(1 



TT 



+ 



^QCD 



rUb 



-2v ■ I 



-2e r 



+ ln^ 



-2v ■ I 



+ 



12 



(16) 



-2e 



In the expression for ^'^^ the first bracket contains the vertex correction for the weak 
current and the contribution from wave-function renormalization for the heavy quark, the 
second term corresponds to the vertex correction for the electromagnetic current and the 
self-energy insertion on the intermediate quark propagator, and the last term corresponds 
to the box diagram. Note that the various terms depend on three different mass scales: 
mb (hard), 2E^l+ ~ mf,A (hard-collinear), and / ~ A (soft). Whereas the first two scales 
are perturbative, the contribution from the box diagram is dominated by soft physics. 
Indeed, power counting shows that the hard and collinear contributions to the box graph 
only contribute at subleading power in A/rrib, as has been observed previously in [lOj. 
Note that the box contribution involves components of the soft spectator momentum 
other than j^]. However, this dependence will cancel in the matching. Finally, it is 
worth emphasizing that there is no leading contribution to the decay amplitude from the 
region of loop momenta that are collinear (not hard-collinear) with the external photon 
momentum. If it were present, this would imply that non-perturbative effects related to 
the hadronic structure of the photon would contribute at leading power, invalidating the 
factorization formula (^. 

In the next step we evaluate the corresponding contributions at one-loop order in 
the effective theory. Generically, they are of the form Ci_ioop ® (OscET)tree + Ctree ® 
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Figure 3: One-loop diagrams in the effective theory whose contribution to the 
amphtude needs to be subtracted in the calculation of the Wilson coefficients. 



(C'scet)i-Ioop- The second contribution, which involves one-loop matrix elements of 
SCET operators convoluted with tree-level coefficient functions, must be subtracted 
from the amplitudes obtained in the full theory. The relevant diagrams are shown in 
Figure 01 We find 



SCET 



1 , 2 f-2v-l\ 777 

- - + InM — + 



^ e2 ■ V U / ' 12 

SCET 



A™=0. (17) 

The first term in the expression for yl^CET corresponds to the first two diagrams in the 
figure, while the second term is obtained from the last graph. Note that this is precisely 
the same contribution as obtained from the box diagram in the full theory. 

The difference of the expressions given in fll6|) and (fT7|) determines the NLO con- 
tributions to the Wilson coefficient functions. Subtracting the pole terms in the MS 
scheme, we find 

Ci = —2 In h (5 — 4 m a;-y) In h m 2 In x.y 

/i ^ fj,-^ 

2 — 3x TT^ 

+ - ^ lna;^-2L2(l-x^) -7- — , (18) 



C2 



Our result for the Wilson coefficient Ci agrees with the expressions for the hard-scattering 
kernel given in pjOjl and JH]; however, eq. ()10j) shows that identifying the kernel with 
Ci misses the large logarithms encountered when the HQET matrix element F{fi) is 
related to the physical S- meson decay constant. (In other words, these papers implicitly 
assume an mfe-dependent definition and normalization of the B-meson LCDA, which 
is unconventional. In our approach the LCDA is normalized to unity, and all mass 
dependence is explicit in the hard-scattering kernel.) Including these corrections and 
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simplifying the answer, we obtain for the hard-scattering kernel at NLO the final result 



T{l+,E^,mb, ^) = 1 + 



_21n2^ + 21n^ + ln2^^ 



2^2(1 -x^) - 5 



77 

T 



(19) 



Note that, despite appearance, the scale dependence of the expression in brackets is 
of the form 21n^(/i/Z+) — 21n(/i//+) + //-independent terms, indicating that it can be 
canceled by the scale dependence of the LCDA under the convolution integral 

in This would be impossible if the argument of the In^-term depended on the hard 
scales or nib. 



4 RG evolution and resummation 

The hard-scattering kernel T contains the logarithms \n{2E^/ n) and ln(2i?^ l^/fi"^), which 
cannot be made small simultaneously for any choice of the renormalization scale /x. 
While for realistic values of the photon energy {E^ ~ mb/2 ~ 2.5 GeV) these logarithms 
are numerically not too large, it is conceptually interesting to gain control over the 
perturbative expansion of the kernel by summing the various logarithms to all orders in 
perturbation theory. 

It was proposed in ^U] to choose the renormalization scale of order /i^ ~ nib A, and 
to identify the remaining large logarithms in (jlSp with those arising in the matching of 
heavy-collinear current operators onto SCET. The precise nature of this identification 
beyond one-loop order was however left unclear. In the spirit of an effective-theory 
approach one would rather prefer to take the renormalization scale further down to a 
value fi = few x Aqcd independent of the 6-quark mass and the photon energy, yet large 
enough for a perturbative treatment. In this way all dependence on the large scales nib 
and E^ becomes explicit and is contained in Wilson coefficients of the effective theory. 

To gain full control over the large logarithms in the perturbative expansion of the 
kernel we perform the matching onto SCET in two steps [HIEICSI- In the first step the 
off-shell fluctuations of the heavy b quark are integrated out by matching onto HQET. 
Hard-collinear modes with momenta scaling like k_, k±) ~ {A, nib, Vi^bA) are inte- 
grated out in a second step. In contrast with .13^, we avoid the explicit construction of 
the intermediate effective theory. Since this theory is needed only for RG improvement, 
it suffices to perform the matching diagrammatically using the method of regions [23| I24j . 
This approach allows us to compute the functions H and J in systematically order by 
order in perturbation theory. Specifically, because the intermediate mass scale of order 
nibA always arises from a scalar product of a soft momentum with a collinear momentum, 
one can perturbatively match onto the intermediate theory by simply setting the soft 
momentum / = 0. In that way only hard fiuctuations associated with couplings to the 
heavy quark are integrated out. This yields precisely the function H. The jet function 
then follows from the ratio J = T/H. 
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Let us illustrate this for the case at hand. For / = the second and third diagrams in 
Figure 121 involve scaleless integrals that vanish in dimensional regularization, while the 
box graph can readily be shown to vanish at leading power. The remaining contribution 
from the weak vertex correction and wave-function renormalization for the heavy quark 
yields 



^QCD 



1=0 



'2t 



1 5 2 In 



2e 



-2 In"' x^-F 



2 3^^-. 



In x^—2L'2{l—x^ 



TT 

12 



(20) 

while the expression for is the same as that for given in (|TH|l . All graphs 

in the effective theory involve tadpole integrals and vanish in dimensional regularization. 
Hence, after MS subtractions the above result determines the hard function Hi defined 
in (fTTj) . Using the first relation in ()13p. it then follows that 



/i 



1 + 



4:71 



-2 In' 



2E^ 



2 In 



7T 



/i 



-2L,il-x,)-A-- 



(21) 

According to the difference between ()19p and (jUJ) determines the one-loop contri- 
bution to the jet function, which is thus given by 



J 



2E^ 



1 + 



47r V 



TT 

"6 



(22) 



In order to proceed we need RG equations obeyed by the various coefficient functions. 
The fact that the decay amplitude in is scale independent links the scale dependence 
of the hard-scattering kernel to the evolution of the LCDA (j)^{u,iJ,) [20]. By analyzing 
the renormalization properties of this function, one finds that the hard-scattering kernel 
satisfies the integro-differential equation (for > 0) [25^^ 



d 



din// 



T(/+,/i) 



rcusp(as) Iny^ + 7(a 
t-i- 



T(/+,/i)+ / dujl+T{uj,l+,as)T{uj, n) 





(23) 

where Fcusp is the universal cusp anomalous dimension familiar from the theory of the 
renormalization of Wilson loops [T^. The function F obeys J dcj F(u;, cj', a^) = 0. At 
one-loop order, it is given by 



t1 — loop ; 

cusp ' 



uj{uj — Uj') Uj'{iO' — Uj) 



(24) 



where the plus distribution is defined such that, when F is integrated with a function 
/(a;), one must replace f^u) — > f{u) — f{uj') under the integral. It is straightforward to 
check that our one-loop result in (fT^ is a solution to at order a^. 

^For simplicity of notation, we omit the arguments E^, mi,, and for the remainder of this section. 
Also, unless otherwise indicated, as = as{ij). 
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From the factorization property of the hard-scattering kernel exhibited in Q and 
the functional forms of the hard and jet functions given above, it follows that the hard 
component and the jet function obey the RG equations 



^ H{^^) 



(iln/i 
d 



-rcusp(a.) In ^ + 7(as) - 7'("s) 



(iln/i 



rcusp(as) In ^^-^ + 7'(tts) 



H{fi) , (25) 

POO 



The reasoning here is analogous to an argument presented by Korchemsky and Sterman 
in their discussion of the B Xs'~f photon spectrum The anomalous dimensions 7 
and 7' do not have a simple geometric interpretation and must be determined by explicit 
calculation. From ()21|) and (j^2I) we find 

7(«.) = -2Cp ^ + 0{aj) , 7'(«s) = 0{a^) . (26) 

From the discussion above, it follows that in the intermediate theory the amplitude 
is represented in terms of the time-ordered product of the electromagnetic current with 
a heavy-collinear SCET current operator of the type x^ih, where x is a hard-coUinear 
quark field. The corresponding matching relation for such currents reads j2] 

u^il - l,)h - Cf'''^{^^) X7"(l - l,)h + . . . , (27) 

where the ellipses represent terms with different Dirac structure. The Wilson coefficient 
Cf^^^{^) coincides with our function Hi{fi). The first relation in (j^^ then determines 
the exact form of the RG equation obeyed by the Wilson coefficients of heavy-collinear 
currents in SCET, which we have derived here for the first time. Based on a comparison 
of one-loop results obtained in SCET with expressions presented in 26^ for the photon 
energy spectrum in inclusive B — >■ Xg'y decays, previous authors have conjectured a 
relation between the coefficient of the ln(/i/2£'^) term in ()25|) and the cusp anomalous 
dimension ^1121; however, the nature and origin of this connection beyond one-loop order 
was left unclear. Note, in particular, that the minus sign in front of the logarithm in 
the RG equation for H{fi) does not allow for a simple interpretation in terms of a cusp 
singularity of a Wilson loop. Our derivation establishes two important facts: First, to 
all orders in perturbation theory only a single logarithm of the ratio fi/2E^ appears 
in the RG equation. Secondly, the coefficient of the logarithm equals minus the cusp 
anomalous dimension. The first observation is a crucial one. The fact that no higher 
powers of logarithms appear in the RG equations cannot be inferred from a fixed-order 
calculation in SCET. It follows from remarkable properties of Wilson loops discussed long 
ago by Korchemsky and Radyushkin p!^. Without this insight it would be impossible 
to integrate the RG equations. The second observation implies that the coefficient of 
the logarithmic term is known to two-loop order, which allows for the resummation of 
Sudakov logarithms at NLO. 
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We now discuss the general solution of the evolution equations ((221) and (j2Sl), which is 
non-trivial due to the convolution integral involving the function T{uj, u', ag). Note that 
the kernel F is not simply a function of the difference {uj — uj') , and so the convolution 
cannot be turned into a product using Fourier transformation. Nevertheless, an exact 
solution can be written down in terms of a new function |2S] 



At one-loop order we find from 



duj uj' r(uj,uj' ,as) ( — 



(28) 



3^ 



1— loop , 



a, a. 



t1— loop f 
cusp 



(a,) ^(1 + a) +^(l-a) +27E 



(29) 



where ip{z) is the logarithmic derivative of the Euler F- function. We start by solving 
the first equation in (f^K|) with the initial condition for H{fih) evaluated at a high scale 
fih ~ ^b, for which it does not contain large logarithms. We then evolve the function 
H{fi) down to an intermediate scale /ij ~ i/mfoA and multiply it by the result J(/+, fii) 
for the jet function, which at the intermediate scale is free of large logarithms and can 
be written in the general form J(/+,/ij) = ^[as(/ij), ln(2_E^ Z+//i^)]. This determines 
the kernel T(/+,/ij) at the intermediate scale. Finally, we solve and compute the 
evolution down to a low-energy scale /i ~ few x Aqcd- The exact solution is given by 



rj=0 



(30) 



where the notation ^[^^(/ij), V,,] means that one must replace each logarithm of the 
ratio 2E^ i+/f^i by a derivative with respect to an auxiliary parameter rj. The evolution 



function U is given by 



as (Mi) 



U{l+,fx,fj,i,(Xh,v) 



r 

, cusp 

da 



» 



/3(a) 



da' 



fJ'h J f3{a') 



da 



da ,^™^P^") 



f3{a) 











f da' 




In^H 


J f3{a') 




/i 




a 





da 



(3{a) 

7(«) 
(3{a) 



2E^ 

+ r]in ^ h 



da ( , f r ,rcusp(a' 



(31) 



where (3{as) = das/d\nfi. Close inspection shows that the result for the hard-scattering 
kernel is independent of the two matching scales and /ij. 

Given this exact result, it is straightforward to derive approximate expressions for 
the kernel at given orders in RG-improved perturbation theory, by using perturbative 
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expansions of the anomalous dimensions and /9-function to the required order. Unfortu- 
nately, controlling terms of 0{as) in the evolution function U would require knowledge 
of the cusp anomalous dimension at three-loop order (as well as knowledge of all other 
anomalous dimensions at two-loop order), which at present is lacking. We will, however, 
control the dependence on the variables 1+ and to 0{as)- As usual, we write 



Tcusp ('3^s 



n=0 



n=0 



(32) 



and similarly for the anomalous dimensions 7 and 7'. The relevant expansion coefficients 
are To = ACp, T, = ACp [(f -^)Ca-^ TpUj] , 70 = -2Cf, 7o = 0, and /3o = ^ 



^TpUf, (3i = C\ — ^ CaTfUj — ACpTprif. Defining the ratios ri = as{^-i) I dsil^h) and 



4 ^ /? _ 34 ^2 20 
f2 = (y.s{ij) / Cis{iii), we obtain our final result 

clnr2 /oT-i \ — clnr 



T(/+,/i) = e 



I, 



X O + 



+ 



4:71 

Air 



f2E. 



, r. 2E^ , 2E^ 
- 2 In^ — ^ + 2 In — ^ 



In- 



2E^ Ia. 



X '■y In X 1^ 



-2Lo{l-x^-A~ 



12 



— + clnr2) — — clnr2) — 27^ 



- ilj'il + clnr2) + ij'il - clnr2) - 1 



TT 



(33) 



where c = Vq/2(3q, and 
To 



£1 _ A 

ro /So 



47r 



t/o(yU,/ii,/i/i) 



n - Inri) 



47r 



+ (l + lnr2) 



47r 



47r 



+ A(inV, + ln^r2)+(^-|) (in !i + 2 - 

7o , . X 1 r(l + clnr2) 

- 77^ ln(rir2) - In — - 2-fE c In r2 + O (a, 

2po r(l-clnr2) 



ri 



r2 



(34) 



corresponds to the evolution function U evaluated with = 0, /+ = /i, and 2Er^ = ^h- 
The only piece missing for a complete resummation at NNLO is the 0{as) contribution 
to Uq, which is independent of /+ and E^. 

In order to study the importance of RG improvement and Sudakov resummation, we 
compare in the left-hand plot in Figure ID the result for the resummed hard-scattering 
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/+ [GeV] /+ [GeV] 

Figure 4: RG-improved predictions for the hard-scattering kernel at maximum 
photon energy. Left: Resuhs at /i = 1 GeV. The bands refer to different values 
of the intermediate matching scale: fif = Ahrrih (center), 2A/jmb (top), O.dAhmi, 
(bottom). Their width reflects the sensitivity to the high-energy matching scale 
l^l, varied between 2ml 0.5ml. The dashed line shows the result obtained 
at one-loop order. Right: Dependence of the kernel on the renormalization 
scale /i, varied between 0.75 GeV and 2.0 GeV as indicated on the curves. 

kernel in ()33|) with the one- loop approximation in ()19|) . We choose the highest possible 
value of the photon energy {E^ = mh/2 with m^ = 4.8 GeV) so as to maximize the 
values of the large logarithms, and plot the function T(/+, /x) for /i = 1 GeV and different 
choices of the matching scales fit and /Xj. Here and below, the "natural" choices fih = 
2E^ and fii = ^j2K^^Kh, where A/j = 0.5 GeV serves as a typical hadronic scale, are 
taken as default values. We use the two-loop running coupling normalized at asimi) = 
0.22 and set = 4 for the number of light quark flavors. (For simplicity, we do 
not match onto a three- flavor theory even for low renormalization scales.) We find 
that resummation effects decrease the magnitude of the radiative corrections, i.e., the 
resummed kernel is closer to the tree- level value (T = 1) than the one- loop result. The 
fact that after Sudakov resummation the radiative corrections are moderate in magnitude 
persists even for asymptotically large 6-quark masses. For instance, setting mi, = 50 GeV 
we find after resummation T(/+, /z) = 0.74 at /+ = = 1 GeV. (Fixed-order perturbation 
theory breaks down for such large values of the quark mass. From (fTIHl we would obtain 
T(/+,/i) = 0.08 with these parameter values.) The figure also exhibits that our results 
are stable under variation of the two matching scales. Varying fif and ^\ by factors of 
2 changes the result for the kernel by less than 10%. This suggests that the unknown 
NNLO corrections to the function Uq in (jH^ are perhaps not very important. 

The scale dependence of the resummed expression for the kernel is illustrated in 
the right-hand plot in Figure EJ which shows the functional dependence of T{l^,fi) for 
maximal photon energy and several values of fi. The matching scales are set to their 
default values fit = = 4.8 GeV and /ij = VA/jm^ ~ 1.55 GeV. We observe a significant 
scale dependence of the kernel, especially as one lowers fi below the intermediate scale 
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[GeV] 

Figure 5: Energy dependence of the convolution integral I{E^) normalized 
to its tree-level value, assuming the model (jHUjl for the LCDA at a low scale 
/io such that as(/io) = 0.5 (top), 0.75 (center), and 1.0 (bottom). The match- 
ing scales are set to their default values. The solid curves correspond to the 
resummed kernel, while the dashed ones are obtained at one-loop order. 



/ij. In other words, the second stage of running (for < /ij), which we have computed 
for the first time in the present paper, is numerically significant. 

Finally, it is interesting to study resummation effects for the convolution integral 

I{E,) = Hdl^ 'MllA T{1^^ E„ m,, /i) . (35) 
Jo 

At tree-level, this integral has been denoted by Jq = 1/Xb, where As is a low-energy 
hadronic parameter [Hj. To evaluate the integral beyond tree-level one needs to assume 
a particular form of the LCDA. In [20], a model function was derived from a QCD sum- 
rule analysis of the matrix element of the bilocal HQET operator in This study 
motivated the ansatz 

7 2 

(pBf^l^^^^) = +e-^+/^B ^ ^ith Xb = -{niB -nib) ^0.32 GeV . (36) 

3 

Such a nonperturbative calculation, which does not control the scale dependence of the 
LCDA, is reasonable only at a low hadronic scale /io- For larger values of n evolution 
effects introduce a radiative tail in the LCDA. We thus expect significant modifications 
of the convolution integral due to radiative corrections. While the function I{E^) is 
formally /i- independent, our model estimate will depend on the value of the scale /io 
at which the functional form given in (j36|) is assumed to be correct. In Figure El we 
show results for the function I{E^) in units of Jq for three different values of yUo- After 
RG resummation we observe a modest reduction of I{E^) with respect to its tree-level 
value, which is fairly insensitive to the precise value of /iq and only shows a mild energy 
dependence. In a rough approximation we have I{E^) ^ 0.75Jo. In contrast, the results 
obtained at one-loop order are strongly sensitive to the choice of /io and exhibit a more 
pronounced dependence on the photon energy. 
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5 Conclusions 



We have applied soft-coUinear effective theory to prove a QCD factorization formula for 
the radiative semileptonic decay B 'yli/, stating that at leading power in A/mb the 
decay amplitude can be written as a convolution of a perturbative hard-scattering kernel 
with the leading-order 5-meson light-cone distribution amplitude. Besides deriving the 
precise form of this convolution and establishing that the kernel is infrared finite to 
all orders in perturbation theory, we have shown that the convolution integral is free of 
endpoint singularities, and that non-valence Fock states of the B meson do not contribute 
at leading power. 

We believe that the analysis of the factorization properties of the decay amplitude 
is most transparent in the framework of the formulation of soft-coUinear effective theory 
developed in 8j, where full QCD is directly matched onto a low-energy effective theory 
in which only long-distance modes are kept as dynamical degrees of freedom. In the 
present case these are the soft constituents of the B meson and the coUinear photon 
field. For more complicated processes with energetic, light final-state mesons one would 
also have to introduce coUinear quark and gluon fields. Many of the techniques used in 
our analysis (such as soft-collinear gauge invariance, reparameterization invariance, and 
renormalization-group improvement) are equally relevant in such a more general context. 

The second part of our analysis was devoted to the calculation of the hard-scattering 
kernel in the factorization formula using renormalization-group improved perturbation 
theory. We have established a second, perturbative factorization formula, according to 
which the different short-distance scales entering in the calculation of the kernel (hard 
scales of order rrib and hard-collinear scales of order y/rribA) can be separated into a hard 
function and a jet function. The corresponding two classes of large logarithms can be 
systematically resummed by solving evolution equations derived from the renormaliza- 
tion properties of the leading-order 5-meson light-cone distribution amplitude. As a 
byproduct, we have elucidated the relation between the anomalous dimension of heavy- 
coUinear currents in the effective theory and the universal cusp anomalous dimension 
encountered in the study of Wilson loops with light-like segments. In contrast to previ- 
ous analyses, we have performed a complete resummation of Sudakov logarithms down 
to a low-energy scale fi ~ few x Aqcd independent of the heavy-quark mass. Only in 
that way all dependence on and the photon energy is explicitly contained in the 
hard-scattering kernel. Our results (j3U|) and (j3H) give the exact analytic solution for the 
kernel, valid to all orders in perturbation theory. 

The discussion of the decay B 7/1/ presented here can be taken over almost ver- 
batim to analyze related processes such as -B ^ 77 and B --fl^l^. The interesting 
observation that at leading power the long-distance effects in these processes are univer- 
sal (for fixed photon energy) ^2] can be understood in our formalism as follows: Once 
short- distance modes related to hard and coUinear interactions with the heavy quark are 
integrated out at a scale /x/i ~ mt, the evolution of the resulting operators in the effective 
theory is independent of the Dirac structure of the relevant current operators. This is 
true for the evolution of heavy-coUinear currents (running from fi^ down to an interme- 
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diate scale /Xj ~ Vm^A), and for the bilocal heavy-light currents whose matrix elements 
are expressed in terms of the leading-order 5-meson light-cone distribution amplitude 
(evolution from /ij down to a low-energy hadronic scale). Process-dependent corrections 
arise only in the initial matching at the high-energy scale and thus are calculable in an 
expansion in powers of as{mb). 

From a phenomenological point of view, an important finding of our analysis is that 
Sudakov resummation does not lead to a strong suppression of the decay amplitudes. 
After renormalization-group improvement we find moderate corrections to the kernel 
which are smaller than at one-loop order, and which stay at the level of 20-30% even in 
the hypothetical case of asymptotically large heavy-quark mass. 
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